A critical reexamination of the non-reciprocal x-ray gyrotropy in V2O3. 
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We critically review the experimental determination of the nonreciprocal x-ray gyrotropy in the 
monoclinic antiferromagnetic insulating phase of V2O3 highlighting some of its paradoxical conse- 
quences. We also note the apparent absence of a normal (non magnetic) linear dichroism in the 
data and speculate on the possible reasons for it. Given the importance of the implications of a 
non reciprocal effect for the physics of V2O3, one of the paradigmatic strongly correlated system, 
we point up the need for at least repeating the experiment under better controlled conditions and 
suggest possible tests for establishing unambiguously its existence. 
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I. INTRODUCTION 



The crystal and electronic structure of V2O3 has been 
the subject of an intensive experimental and theoretical 
study throughout the 70's 1-6 that led to the identifica- 
tion of V2O3 as the prototype of Mott-Hubbard systems. 
This compound shows a metal-insulator transitions, due 
to the interplay between band formation and electron 
Coulomb correlation, from a paramagnetic metallic (PM) 
phase at room temperature to a paramagnetic insulating 
(PI) phase at higher temperature (w 500 K) and to an 
antiferromagnetic insulating (AFI) phase at lower tem- 
perature (T c « 150 K for the stoichiometric compound). 
Associated with this latter there is also a structural tran- 
sition from a corundum cell (rhombohedral system) to a 
monoclinic one (see Fig. 1). 

In this investigation a correct description of the prop- 
erties of the AFI ground state has been of paramount im- 
portance for both the physical understanding of the insu- 
lating phase and for throwing light on the mechanism of 
the metal insulator transition. Indeed Castellani, Natoli 
and Ranninger (CNR), 5 ' were the first to realize the im- 
portance of the orbital degrees of freedom for the expla- 
nation of the peculiar antiferromagnetic spin structure 3 
that breaks the trigonal symmetry of the corundum lat- 
tice in the high temperature phase. Evidence for their 
role was subsequently found in a series of inelastic neu- 
tron scattering experiments by Wei Bao et al. 7 and NMR 
nuclear spin relaxation measurements 8 in the high tem- 
perature paramagnetic phases. The ensuing physical pic- 
ture was one in which the orbital degrees of freedom are 
frozen in the AFI phase, giving rise to the peculiar spin 
structure, and are responsible for the short range mag- 
netic fluctuations in the high temperature phases, provid- 
ing much of the necessary entropy content for the transi- 
tion to occur. Even though the Hubbard model used by 
CNR was based on a spin S = 1/2 for the V 3+ ion in a 
doubly degenarate orbital e g state, the subsequent pic- 
ture based on spin S=l, implied by some experiments 9,10 
and worked out by Mila et al. 11 ' 12 and Di Matteo et al., 13 



did not change the above physical description, but only 
redefined the basic building blocks of the problem. In- 
deed in the S=l picture the elementary units are the ver- 
tical biatomic molecules in a doubly degenerate ground 
state, strongly coupled to spin S=2, and weakly interact- 
ing via electronic kinetic hopping in the basal planes. 

The introduction of the extra degrees of freedom in the 
problem might be expected to change the magnetic space 
group of the ordered AFI phase as inferred from the lat- 
tice geometry and spin moments only. As explained in 
the following section, the magnetic space group of the 
monoclinic cell with the atomic sites occupied by the 
magnetic moments in the way found by Moon 3 is the 
magnetic space group P2/a + T{E\to}P2/a, 14 contain- 
ing explicitly the time reversal operator T and the inver- 
sion I (E is the identity operator, to is the body-centered 
translation) . 

In this context the observation of a nonreciprocal x-ray 
linear dichroism measured at the K-edge of the Vanadium 
atom in the AFI phase of V2O3 15 seemed to shed new 
light on the properties of the insulating ground state, 
since the effect can only be observed when neither the 
inversion nor the time reversal are separately symmetry 
operations of the system, thus implying that the system 
should be magnetoelectric. Even though not explicitly 
stated in Ref. [ 15] this finding was pointing to a role of 
the orbital degrees of freedom in lowering the magnetic 
symmetry of V2O3, all the more that an independent 
observation of the forbidden (lll) m reflection in the V 
K-edge anomalous diffraction 9 had been interpreted as 
implying some kind of orbital ordering. 

However a closer consideration of all the pieces of evi- 
dence of the problem revealed some inconsistencies in the 
overall interpretative frame, as discussed at length in Ref. 
[13]. First of all the theoretical model used in Ref. [ 9] to 
interpret the orbital ordering was based on the old CNR 
S=l/2 model, which was later shown to be inadequate 
to describe the physics of V2O3. In the new correlated 
model based on spin S=l 11-13 the orbital ordering in the 
AFI ground state was again found compatible with the 
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P2/a + T{E\to}P2/a group. In the meantime growing 
evidence was accumulating that the forbidden (lll) m re- 
flection is of pure magnetic origin and has nothing to do 
with orbital ordering. 16-18 Last but not least, a new di- 
rect measurement of the magnetoelectric effect in V2O3 
has provided a negative result, 19 confirming the previous 
result by Astrov. 20 

Concomitant to all this, our efforts to understand at a 
deeper level the implications of the nonreciprocal x-ray 
gyrotropy experiment led us to note some paradoxical 
consequences of this result and to look for alternative in- 
terpretations. It is the purpose of this paper to present 
the conclusions of this investigation and to point up at 
least the need to repeat the experiment under better con- 
trolled conditions. 

This paper is organized as follows. In section II, we dis- 
cuss in some detail, for the convenience of the reader, the 
structural and magnetic properties of the system in its 
monoclinic phase and derive the corresponding magnetic 
space group, listing all the magnetoelectric subgroups. In 
section III we derive explicitly the response functions cor- 
responding to a linear dichroism experiment both in the 
dipole-dipole (El-El) and the dipole-quadrupole (E1-E2) 
channel. We show that the corresponding signal is pro- 
portional to the average of certain operators in the final 
scattering state of the excited photoelectron, having def- 
inite transformation properties under inversion and time 
reversal, and discuss in the light of the possible magnetic 
space groups the conditions under which such a signal is 
different from zero. We find that contrary to what im- 
plied in Ref. [ 15] one expects a significant signal from the 
El-El channel, which is inversion and time reversal even. 
We further argue that its existence is a necessary conse- 
quence of the Templeton scattering at the (101) m forbid- 
den reflection in the a — a channel observed by Paolasini 
et al. 21 in the monoclinic phase of the same compound. 22 
We estimate the dichroic signal in the frame of multiple 
scattering (MS) theory, which should be adequate to de- 
scribe transitions to the itinerant conduction p states in 
V2O3. Surprisingly enough its shape is quite similar to 
the "nonreciprocal" gyrotropy signal reported in Ref. [ 
15]. We speculate on such similarity and give arguments 
whereby the reversal of the signal under magnetic field 15 
could in fact be due to an uncontrolled geometrical effect. 
In section IV we summarize our conclusions. 



II. THE MONOCLINIC STRUCTURE OF V2O3. 

Pure V2O3, as well as the 2.8% chromium doped sam- 
ple used in the x-ray experiments mentioned in the in- 
troduction, is a paramagnet at room temperature that 
crystallize in the trigonal system. At T c ~ 150 K (180 K 
for the doped compound), it undergoes a structural tran- 
sition to the monoclinic system and an antifcrromagnetic 
order sets in. This doubles the unit cell to four formula 
units, instead of the two formula units of the trigonal 



cell. Moreover the transition breaks the trigonal symme- 
try in such a way that one of the three originally equiv- 
alent vanadium-vanadium bonds in the hexagonal plane 
of the corundum cell becomes ferromagnetic and longer 
of about 4% than the other two, which are antiferromag- 
netic. Because of this, three possible monoclinic twin do- 
mains exist, depending on which of the three equivalent 
bonds stretches and becomes ferromagnetic. They are re- 
lated by a rotation of 120° around the corundum c^-axis. 
Experimentally one finds that V2O3 has the tendency to 
crystallize prevalently in one of the three monodomains, 
more than in a statistical average of the three and that 
the choice among them is more or less random, depend- 
ing on the size, shape and boundary conditions of the 
crystal. 23 

We give here a detailed description of the monoclinic 
cells and of its lattice and magnetic symmetries, as we 
shall use them in the subsequent analysis of the structure 
factor. 




~-O^Cr~ and have "opposite" oxygen octahedra 

FIG. 1. Corundum and monoclinic cells for the paramag- 
netic and the antiferromagnetic phases of V2O3, respectively. 
The dfr-vector coincides with b m and bn connects 4' with 4. 
The distortion and increase in volume due to the monoclinic 
transition are not shown. Note that V-ions do not lay on 
the hexagonal planes, but they are slightly displaced (0.17A) 
towards the nearest voids. 

The crystal positions are taken from Dernier and 
Marezio, 2 and the magnetic structure of the monoclinic 
cell from Moon 3 later confirmed by Wei Bao. 7 

The monoclinic cell, shown in Fig. 1, belongs to the 
body-centered crystallographic space group I2/a. Us- 
ing the reference frame and the numbering of vana- 
dium atoms of Fig. 1, we can divide the eight atoms 
of the monoclinic cell into two groups of four, with 
opposite orientation of the magnetic moment, V\ = 
(1/2 - u,v,-w), V 2 = (1/2 + u,-v,w), V 3 = 
(—u,—v,—w), V4 = (u,w,w), and Vy, Vy, V3/, V4', with 
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coordinates obtained by adding the vector (1/2, 1/2, 1/2) 
to the first group. Here u = 0.3438, v = 0.0008, w = 

0. 2991 are the fractional coordinates of the atoms in unit 
of the monoclinic axis. 2 Note that, neglecting the mag- 
netic moments, the two groups of atoms with their oxy- 
gen environments are translationally equivalent. 

From these data we can infer that the magnetic space 
group of the monoclinic cell is the magnetic space group 
P2/a + f {E\t a }P2/a, 14 where f is the time-reversal op- 
erator and the monoclinic group P2/a contains the iden- 
tity E, the inversion I, the two-fold rotation about the 
monoclinic b m -axis C' 2b and the reflection rh b with re- 
spect to the plane perpendicular to this axis. Of course, 
the appropriate translation is associated to each of these 
operators. 13 Choosing the origin of the system as in Fig. 

1, we have the following symmetry operations: 



1) e, i 

2) C 2b , m b 

3) f , fl 

4) TC 2bl trhb 



No translation 

2 (Pm ~l~ ) 
1 tz 
2 



(& m -\- b m -\- c m ) 



7^7 1 



Notice that, if one considers only the lattice and mag- 
netic moments, the presence of the time-reversal symme- 
try is a necessary consequence of the translational equiv- 
alence of the two groups of atoms. Now, as anticipated in 
the introduction, one of the consequences of the observa- 
tion of non-reciprocal effects by linear dichroism 15 is the 
indication that other degrees of freedom may play a role 
in determining the symmetry of the insulating ground 
state. In fact a necessary condition to detect the non- 
reciprocal gyrotropy tensor is that the system be mag- 
netoelectric, this fact implying in turn the breakdown of 
time-reversal and inversion symmetry. Orbital degrees of 
freedom are the natural candidates in the physics of V2O3 
to obtain such a symmetry reduction, although new, 
still unsuspected, ingredients might come into play. One 
other possibility, already debated in the literature, 3 ' 24 ' 25 
yet very controversial, concerns the presence of an out-of- 
plane (a m -c m plane) component of the magnetic moment 
directed along the y axis.. This fact would at least break 
the glide-plane symmetry. Three independent neutron 
scattering experiments 3 ' 24 ' 25 were not able either to dis- 
card or accept this hypotesis. However this out of plane 
component would contradict the argument given in Refs. 
[ 13, 16-18] to explain the absence of any dipole signal at 
the (lll) m reflection in the energy region of the conduc- 
tion band of p symmetry, that can only be explained by 
the presence of the glide-plane symmetry together with 
the fact that < L y >= 0. For this reason, we can con- 
clude that x-ray measurements point toward the absence 
of an out-of-plane component of the magnetic moment. 

Here we list for completeness all the possible space 
magnetic subgroups of P2/a + T{E\to}P2/a that do ad- 
mit magnetoelectricity: 



1) {E,fl,fC 2b ,m b ) =P2'/a 

2) {E,fl,C 2b ,fm b ) =P2/a' 

3) (E,fl) 

4) (E,TC 2b ) 

5) (E, C 2b ) 

6) (E,fm b ) 

7) (E,m b ) 

8) (E) 

We have indicated in parenthesis the symmetry opera- 
tions of the associated point groups and for simplicity we 
have not written the corresponding translations which 
are the same as above. Note that all the two elements 
subgroups actually belong to the triclinic system. On the 
basis of these groups we shall discuss in the next section 
the conditions for the observation of linear dichroism in 
absorption. 



III. NORMAL AND NONRECIPROCAL LINEAR 
DICHROISM 

Linear dichroism in the x-ray range is the differen- 
tial spectrum obtained by subtracting two core absorp- 
tion spectra taken with orthogonal linear polarizations. 
According to whether it is of dipole-dipole (El-El) or 
dipole-quadrupolc (E1-E2) origin we shall distinguish, for 
convenience and reasons to become apparent later, be- 
tween respectively: Normal (NXLD) and Nonreciprocal 
(NRXLD) X-ray Linear Dichroism. The two contribu- 
tions can obviously be both present in the same spec- 
trum. The experiment 15 we are interested in was per- 
formed with the x-ray beam directed along the hexagonal 
c/j-axis, or (202) m in monoclinic Miller indices, with the 
two orthogonal directions of the linear polarizations ly- 
ing in the hexagonal plane without specification of their 
orientation with respect to the in plane crystallographc 
axes. 



A. Dichroic operators and related sum rules 

As is well known, at the Vanadium K-edge the x-ray 
absorption cross section is given by: 

at = 4tt 2 ^ £ |<*W \d\^)\ 2 6(hw - (E n - E )) (1) 



The operator O = e ■ f(l + %k ■ r) is the usual 
matter-radiation interaction operator expanded up to the 
quadrupolar term, with the usual notation for the photon 
polarisation e and the wave vector k. is the 

ground (excited) state of the crystal, E (E n ) its energy 
and the index i indicates the lattice site of the Vanadium 
photoabsorbing atom. The sum is extended over all the 
excited states of the system and Tito is the energy of the 
incoming photon, while a is the fine-structure constant. 
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Before analyzing Eq. (1) on the basis of the crystal 
space magnetic symmetry, we derive explicitly the cross 
section and the consequent sum rule for non-reciprocal 
and normal linear dichroism, on the basis of the formal- 
ism developed in Ref. [ 26]. In both cases linear dichro- 
ism is defined as: 



NRXLD 



4ir 2 ahLoM?(E) J dk e 



a i = Gi{c a ) - ai(e b ) 



(2) 



where e a and e b are ususally chosen orthogonal. With 
reference to Sec. 7 of Ref. [ 26], we can write the signal 
up to the E1-E2 channel as the scalar product of two ten- 
sors, representing respectively the properties of the light 
(T) and of the matter (M): 



rt D cx^(-) 1 +^[rf (ea , ea ,fc)M^ 



-T^(e b ,e b ,k)M^} 



(3) 



where p=l,2,3. The explicit expression of the tensors 
Tq P ^ is given in Sec. 7 of Ref. [ 26]. Moreover, in terms 
of the tensors expressing the properties of the matter (M, 
for El-El and M for E1-E2) we obtain: 



NXLD 



NRXLD 



(3) 



oc((M, 

T^ 1 



M 



(3) 



M 



(2) 



C.C.) + 
C.C.)) 



(4) 



(5) 



In both sets of equations we have omitted for simplic- 
ity the index i in the right hand side. To evaluate the 
matter tensor we shall adopt the one particle approach 
in the framework of multiple scattering theory. However, 
as shown in Ref. [ 26], all the expressions remain valid 
in a second quantization scheme, if one takes the multi- 
ple scattering functions as the one particle basis for the 
matter operators. Around each site the photoelectron 
wave function can be expanded in spherical harmonics 
and spin states as: 

|*i(fc e ) >= E B LAke)RLa(r;E)Y lm (f) X , (6) 

where k e is the photoelectron wave vector associated to 
the final scattering state with kinetic energy E oc k 2 and 
Xo is the spin state. Using this expression it is possible 
to write the tensors M q in terms of the scattering am- 
plitudes B\ m(J {k). In the case of interest (K-edge) I is 
either 1 or 2 so that we find in the geometrical setting of 
the experiment of Ref. [ 15]: 



NXLD 



4ir 2 ahwM 2 (E) J dk e ^\ m \( B ima)*Bi-m, 

OC J dk e < *i{ke)\{Ll, - L 2 y ,)\*Se) > (7) 



T,m,a\ m \( B 2ma)* B l, 



C.C. 



oc /<4 < V l (k e )\(L 2 x , -Ll,)Q z \^i{k e ) > -c.c. (8) 

where, indicating by R mt the muffin-tin radius, the pro- 
portionality factor is 



M t (E) 



J H ^r*drR?(r;E) 



In the previous equations we have ignored the weak 
ma dependence of the radial wave functions, as justified 
for deep core (Is) electron transitions. We have intro- 
duced the dipole radial transition matrix element, Mi(E) 
and integrated over all the direction of the escaping pho- 
toelectron (k e ) at a fixed photon energy, as appropri- 
ate for absorption. 26 Moreover we have used the equality 
(L 2 — L 2 ) = (L\ + L 2 _) and the fact that the magnetic 
quantum number m runs from -1 to 1. The labels x' and 
y' indicate two arbitrary orthogonal directions of the two 
polarizations e a and e b in the hexagonal plane. As dis- 
cussed in Ref. [ 26] the energy integrated spectrum pro- 
vides the expectation value of the same operators over 
the ground state, a situation that is in common with all 
kinds of dichroisms. 

Equations (7) and (8) show that linear dichroism at 
the energy Uui of the incoming photon measures the 
expectation value of a physical operator over the final 
scattering state corresponding to that energy. In par- 
ticular, this operator is proportional to L 2 , — L 2 , in 
the dipole-dipole channel and to (L 2 , — L 2 ,)Q Z — c.c. in 
the dipole-quadrupolc channel. In the last expression 
fl z = (Lx n — nx L) z , where L is the usual angular mo- 
mentum operator and h the radial unit vector operator. 
It is the raising/lowering operator for the orbital quan- 
tum number I when applied to a spherical harmonic with 
the same m projection, 27 and is known as the anapole 28 
or toroidal moment operator. 

Note that the linear dichroism signal is time-reversal 
even and inversion even in the dipole-dipole channel, 
while it is time-reversal odd and inversion odd in the 
dipolc-quadrupole channel, as expected on the basis of 
general considerations. 



B. Extinction rules for linear dichroism. 



In the monoclinic crystal, the total absorption signal 
is obtained by summin| over all the eight inequivalent 
vanadium sites Ic — Y]j—i 0» m the unit cell. Therefore, 
indicating with the label LD both types of dichroism we 
obtain 
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lLD = j2°f D (9) 

i=l 

By exploiting the symmetry operations of the system, 
we can now relate one another the various af D in the unit 
cell, exactly as done for the anomalous scattering ampli- 
tude in Ref. [13]. Absorption is in fact the imaginary 
part of the forward scattering amplitude. A symmetry 
operation then acts both on the site index and on the 
operator appearing in the averages (7) and (8), assuming 
that the final scattering state |*j(fc e ) > does not break 
the symmetry of the ground state. This latter assump- 
tion is physically plausible since in the one electron ap- 
proximation the rest of the system is seen by the excited 
photoelectron as a static scattering potential. Note that 
the symmetry considerations in this subsection do not 
depend on the one particle approximation introduced in 
the previous subsection, but hold a more general validity, 
relying only on the symmetry group of V2O3.. 

The site transformation rules for the largest magnetic 
group P2/a + T{E\t }P2/a are given in the following 
table: 



E : 


01 


02 


03 


04 


01' 


02' 


03/ 


4 ' 


I : 


C2 


01 


04 


0"3 


02' 


01' 


04' 


03' 


Cib ■ 


(T 4 / 


03' 


02' 


01' 


04 


03 


02 


01 


m b : 


3 , 


cr 4 ' 


01' 


2 ' 


03 


4 


01 


02 


T : 


ay 


0- 2 ' 


3 ' 


0- 4 ' 


01 


02 


03 


04 


TI : 


02' 


01' 


04' 


03' 


C2 


01 


04 


03 


TC 2b : 


CT 4 


03 


02 


01 


04' 


03' 


02' 


01' 


Tm b : 


03 


04 


01 


02 


03/ 


4 ' 


01' 


02' 



Therefore, eg, Taf RXLD = -<r y JRXLD and similarly for 
all the other cases. 

Using these relations we can then express I^d in term 
of one or two absorption site according to the symmetry 
group chosen. We consider the three cases of the groups: 
P2/a + f{E\t }P2/a, P2'/a and P2/a' which are perti- 
nent to the monoclinic phase. 

In the case of the original group P2/a + T{E\t }P2/a 
(non ME), the global signal can be expressed in terms of 
only one independent center of absorption. Thus Eq. (9) 
can be written as: 

Ild = (l + fi)(l+fm b )(l + f)a% D (10) 

In the case of a ME subgroups we use only the part of 
the table referring to the operators of each subgroup so 
that 

I LD {P2'/a) = (1 + TJ)(1 + fm b )(^ D + 4 L , D ) (11) 
and 

I LD {P2/a') = (1 + T7)(l + m b )(a% D + 4 L , D ) (12) 

in terms of two independent centers of absorption. In 
these formulas it is now intended that the group opera- 
tors have already acted on the site but not on the physical 
observables. 



In all cases it is straightforward to derive that no cir- 
cular dichroism can be detected for the system, as exper- 
imentally found in Ref. [ 15], since all the expressions are 
proportional to (1 + T7) and the product of time-reversal 
and inversion symmetries is always -1 for circular dichro- 
ism (1 = 1 and T = — 1 in the El-El channel and I = — 1 
and T = 1 in the E1-E2 channel). As expected, the 
three magnetic groups present different extinction rules 
for NRXLD. Since in this case the relevant operator is 
time-reversal and inversion odd, the signal is zero for the 
non ME group Eq. (10), due to the term (1 + T), and 
different from zero for the other two, since TI = +1. 

On the contrary there is no case in which the linear 
dichroism in the El-El channel (NXLD) is forbidden by 
the crystal structure, because in this channel both the ac- 
tion of the time-reversal and the action of the inversion 
operators give +1. Thus, since the magnetic symmetry 
around each vanadium ion is C\ (ie, no symmetry at all) 
a non zero atomic NXLD is in principle expected. Such 
a dichroic signal is zero only when each pair of vanadium 
ions contribute to the total intensity with an opposite 
sign. This latter occurrence is possible only for the par- 
ticular direction at an angle 7r/4 with respect to the glide 
plane, in which case for each pair of vanadium ions re- 
lated by this symmetry operation (eg 1 and 3', 2 and 4', 
3 and 1', 4 and 2', as inferred by the transformation ta- 
ble given above) L 2 X — > L 2 y and viceversa, leading to the 
extinction of the total crystal signal as seen from Eqs. 
(10), (11) and (12). Note that any reduction of the sym- 
metry to a group that does not contain the glide plane 
(cither multiplied by the time-reversal or not) implies 
that the NXLD signal is never zero. Note also that the 
non-reciprocal signal, proportional to (L\ — Ly)Q z — h.c, 
at the particular angle tt/4 vanishes if the magnetic group 
contains the glide-plane alone as a symmetry operation 
(like P2/a!) and may be different from zero with groups 
non containing it (like P2'/a). 

We dwelt a bit on this analysis, since it highlights an 
unexpected feature in the nonreciprocal XLD spectra as 
presented in Ref. [15]. In fact, since one cannot dis- 
entangle the El-El from the E1-E2 channel one would 
expect the superposition of two linear dichroic signals in 
any experimental spectrum, one of which is time-reversal 
and inversion even (coming from the El-El transition) 
and a second, which is time-reversal and inversion odd 
(coming from the E1-E2 transition). As a consequence 
the residual signal obtained by taking the half sum of 
the two nonreciprocal XLD spectra, respectively with the 
external magnetic field oriented parallcl/antiparallel to 
the hexagonal c H axis, (ie [XLD(H+) + XLD(H")]/2), 
should give the normal (time-reversal and inversion even) 
linear dichroism signal, as correctly recognized by the 
authors in Ref. [ 15]. From Fig. 2 of the same refer- 
ence one can infer that this residual signal is about ten 
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times smaller than the nonreciprocal signal (which is it- 
self between 0.0-1.0 % of the main absorption), in prac- 
tice nearly coinciding with the background noise, 'due 
perhaps to the slightly distorted monoclinic structure', 
in their own words. 

We show in the following section that a multiple scat- 
tering realistic: simulation of the NXLD in the monoclinic 
lattice provides a signal which can be up to 4% of the 
main absorption. The presence of such a signal is also 
corroborated by the sizable Templeton scattering at the 
(101) m forbidden reflection in the a — a channel observed 
by Paolasini et al 21 in the monoclinic phase. 22 We per- 
formed a numerical simulation of the (101) m reflection 
(see next subsection) and found that the signal is almost 
entirely due to the El-El channel. Its existence implies a 
different value of the the x and y dipole matrix elements, 
which is at the basis of the NXLD. Notice that the Tem- 
pleton scattering is zero in the corundum (paramagnetic) 
phase, as experimentally verified in the same work, 21 due 
to the existence of the trigonal symmetry axis C3, since 
then the x and y components of the electron scattering 
wavefunctions are degenerate in the real basis representa- 
tion of the spherical harmonics. As a consequence mixed 
matrix elements of the kind < \I>o|#|*n >< ^nl 2/1*0 >, 
which are at the origin of the Templeton effect in the 
a — a channel at the (101) m reflection (see next subsec- 
tion), vanish. 

The only possibility to have a situation like the one 
depicted in Rcf. [ 15] is in the case that by mere accident 
the two orthogonal linear polarizations of the incoming 
photon beam in the geometrical setting of the experiment 
lay at an angle of 7r/4 relative to the monoclinic b axis 
normal to the glide plane (whether or not multiplied by 
the time-reversal operator). At the same time to avoid 
the extinction of the NRXLD signal, the magnetoclcctric 
group should be P2'/a, since in this case the operator 
Trhb does not change the sign of (L 2 — L 2 )Q Z — h.c. 

Under these assumptions the search for this particular 
geometrical setting would have directly determined the 
ground state symmetry of the AFI phase. However there 
is no trace in Ref. [ 15] of such investigation, so that 
their result cannot be unambiguously interpreted. 



C. Multiple Scattering simulations of NXLD 

The MS simulations of the normal linear dichroism in 
V2O3 were carried out in the framework of the multi- 
ple scattering theory within the muffin-tin approxima- 
tion. We chose a cluster containing 135 atoms, ie 54 
Vanadium and 81 oxygen atoms in the correct ratio 2 
to 3, having a radius of 6.9 A, enough to get conver- 
gence with the cluster size. This latter is assured by 
convoluting the spectrum calculated with the real part 
of the Hedin-Lundqvist (HL) potential with a lorentian 
function having an energy dependent damping T(E) de- 
rived from the universal mean free path curve by the 



relation X(E) = l/k e E/T(E), where E is the photo- 
electron kinetic energy (in Rydbergs) and k e — \f~E~ its 
wave vector. This procedure seems to give better results 
than calculating the spectrum directly with the complex 
HL potential, as recently found in other cases. 29 The ac- 
tual calculations were performed independently with two 
different codes, the FDMNES program 30 which is non 
muffin-tin but incorporates a muffin-tin version and the 
CONTINUUM program, 31 with similar results. As antic- 
ipated in the introduction, the independent particle ap- 
proximation should be adequate to describe transitions 
to dclocalizcd conduction states with p symmetry around 
the photoaborber. Only strong quadrupole transitions to 
empty correlated d states might not be well described, 
together with p transitions to hybridized p — d states in 
the same energy range. However we are only interested 
in the higher energy region where the nonreciprocal gy- 
rotropic effect seems to be more pronounced. Moreover 
in the same energy range the muffin-tin approximation is 
reasonably accurate, due to the close-packed structure of 
V2O3 (despite the Vanadium voids) and the moderately 
high kinetic energy of the excited photoelcctron. In order 
to eliminate any doubt on the confidence of the muffin- 
tin calculations, we also performed non- muffin-tin calcu- 
lations on a substantially smaller cluster (eleven atoms) 
with the FDMNES program 30 and compared with the 
muffin-tin counterpart. General shapes and, more im- 
portant, amplitudes of the various signals discussed be- 
low did indeed compare very favorably. (Larger cluster 
calculations with no local symmetry, like in the present 
case, are probitively long in cpu time and require a lot of 
computer memory). 

Even though the formalism of spherical tensors is very 
convenient to write down the response functions of x-ray 
spectroscopies, as illustrated in section III B, working in 
cartesian coordinates leads to immediately readable for- 
mulas. In the geometry of Ref. [ 15] the direction of the 
incoming x-ray beam coincides with the trigonal axis of 
the high-temperature phase of V2O3. Thus, in the refer- 
ence frame of Fig. 1, the two polarizations have compo- 
nents: e a = (cos a, sin a, 0) and e& = (— sin a, cos a, 0). 

The atomic linear dichroism in the El-El channel is 
then given by: 

a NXLD K 5(u; n ) ((%\e a ■ r1M>„)(*„|e • r\%) 

n 

(%\e b -^ n )(* n \e b -r\%)) 
= ^2S(uj n ) {cos2a((x\x) l - (y\y)i) 

n 

+ sm2aRe[(x\y) l }} (13) 

writing for brevity oj n = fiuo — (E n — Eq), 
(* |x|*„)(*„|a;|* ) — > (x\x)i and similarly for the other 
components. 

When summing over all sites of the unit cell, only the 
part proportional to cos 2a survives, since the mixed ma- 
trix elements Re[(x\y)i] give opposite contributions for 



G 



pairs related by the glide plane. Therefore the crystal 
absorption is given by: 

ag XLD <x cos2a]T<5K)((x|x) - (y\y)) (14) 

n 

This is in contrast to what happens for Tcmplcton scat- 
tering in the a — a channel at the (101) m reflection in the 
same El-El channel. In this case the incoming and out- 
going polarization have components e a = (cos (3, sin (3, 0) , 
so that the scattering amplitude for the atom at site i 
reads: 

fi((j) = {huj) 2 U J 2 ( cos2 @( x \ x )i + sin2 0(y\y)i 

n n 

+ sm2/3Re[(x\y) l }) 

Since the crystal structure factor at the (101) m re- 
flection substracts the amplitudes coming from pairs of 
atoms related by the glide plane symmetry, the net crys- 
tal amplitude Fc(to) turns out to be: 

F C Hcx(M 2 E^ 1 S sin2 ^ e [( a; iy)] ( 15 ) 

n 11 

as anticipated above. Therefore in C3 symmetry both 
linear dichroism and Templeton scattering at the (101) m 
reflection are zero, whereas if this latter does not vanish 
a normal linear dichroism is also expected. 
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FIG. 2. Calculated Templeton intensity at the (101) m re- 
flection. The ordinate shows the square of effective number 
of scattering electrons. 

However it is a bit hard to derive the size of the lin- 
ear dichroism from the size of Templeton scattering, due 
to the extinction corrections in this latter spectroscopy. 
Our aim here is to show that the energy shape of the 
Templeton intensity is mainly of dipolar origin. To this 
purpose we present in Fig. 2 the El-El contribution of 



our simulation for the a — a channel at the (101) TO re- 
flection. The agreement with the experimental spectrum 
of Ref. [21] is rather good, in the sense that we obtain 
peaks and valleys at the right energy position. Only the 
intensities disagree with the experiment. However one 
should observe that fig. (2b) of Ref. [ 21] represents the 
raw data not corrected for absorption. The amplitude 
of the E1-E2 contribution has been calculated to be ten 
times smaller than that in the El-El channel, due to ra- 
dial matrix elements effects and the very small weight of 
the nd (n > 3) scattering amplitudes in the conduction 
band of p symmetry. Since in charge scattering the two 
contributions do not interfere, if one neglects magnetic 
order, the ratio of intensities in the two channels is 1 to 
one hundred. 
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FIG. 3. Absorption signals (Mbarn) in V2O3 correspond- 
ing to polarizations within the hexagonal plane, parallel and 
perpendicular to b m . 

The main results of the numerical calculation for ab- 
sorption and linear dichroism are shown in Figs. 3 and 
4. Figure 3 shows the absorption signals in Mbarn 
with the incoming photon polarizations in the hexago- 
nal plane and directed along b m and perpendicular to 
it. The general shape follows quite well the experimental 
counterpart. 15 In Fig. 4 instead we show on an expanded 
scale the variation of the dichroic signal (in Mbarn) as a 
function of the angle a. The full line curve with negative 
maximum corresponds to a = 0° (< x\x > — < y\y >), 
while the long-dashed line with positive maximum is for 
a = 90° (< y\y > — < x\x >), in keeping with the cos 2a 
variation. The other two curves are for a = 15° and 
a = 75° for later use. It is immedately evident that the 
dichroic signal vanish and then reverts its sign when the 
polarization angle crosses 45°. Its shape is fairly simi- 
lar to the non reciprocal experimental signal reported in 
Ref. [ 15] while the maximum changes from to 4% of 
the main absorption and is about 3% for a = 75°, to be 
compared with the experimental 1%. In both figures we 
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have not corrected for the unknown experimental resolu- 
tion. 

This finding is to be contrasted with an intensity more 
than 10 times smaller for the E1-E2 signal. This latter is 
in fact zero in the paramagnetic phase and picks up inten- 
sity only because of the magnetic ordering. Therefore on 
top of the factor one tenth in the matrix element, men- 
tioned above when discussing Templeton scattering, one 
should add a further reduction factor originating from 
magnetic effects. We did not attempt to calculate the 
atomic signal as given by Eq. (8) for comparison, since 
the global signal depends on the crystal magnetic space 
group, as discussed above. Not knowing the mechanism 
that determines the lowering of the magnetic symmetry 
makes it impossible a meaningful estimate of the crystal 
magnetic effect. However this discussion should indicate 
that the non reciprocal gyrotropic effect in V2O3, if it 
exists, is more than one order of magnitude smaller than 
normal linear dichroism. 
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FIG. 4. Linear dichroism signals (Mbarn) as a function of 
the azimuthal angle a about the wave vector parallel to z. 
The four curves correspond to a — 0, 15, 75, 90° degrees . The 
labeling of peaks corresponds to that of Ref. [15]. 

The relevance of these calculations to the argument we 
are discussing stems from the fact that there is a remote 
chance that the reversal of the dichroic signal with the 
magnetic field might be a geometrical artifact due to a 
lack of full control of the experimental conditions. In 
fact, since the procedure of magnetoelectric annealing is 
rather ineffective due to the non magnetoelectricity of 
V2O3 and the sizable conductivity of the sample in the 
paramagnateic insulating phase as discussed in Ref. 19, 
it is not hard to conceive that the procedure performed 
in Ref. 15 could in reality involve two twins rotated by 
120° degrees, so that the two curves refer to linear po- 
larizations at 15° and 75° degrees with respect to the 
monoclinic b m -axis. Therefore the control of the orien- 
tation of the crystal domain after annealing is essential 



in order to eliminate any ambiguity in the interpretation 
of the results. 

Admittedly this is a quite remote possibility, however 
not more remote than the one depicted at the end of 
the previous subsection under the assumption that the 
observed effect was a truly nonrcciprocal effect. Clearly 
the whole argument hinges on the existence of a sizable 
linear dichroic effect of El-El character, which is anyway 
a prediction of the present work subject to experimental 
test. In fact a measurement of the angular variation of 
the linear dichroism in the AFI phase without magneto- 
electric annealing might settle the question. 



IV. CONCLUSIONS 

In the preceding sections we have argued about the 
existence of a measurable normal [it time-reversal and 
inversion even) linear dichroism in the AFI monoclinic 
phase of V2O3, both on the basis of realistic simulations 
in the framework of MS theory and experimental evi- 
dence coming from the observation of Templeton scat- 
tering at the (101) m forbidden reflection in the a — a 
channel. 21 Therefore it comes to a surprise that no such 
a signal is present in the nonreciprocal linear dicroism 
experimental spectrum measured by Goulon et al 15 in 
the same sample. 22 From the analysis of the crystal ab- 
sorption cross section according to various ground state 
magnetic symmetry groups we have concluded that there 
exists only one group (P2 f /a) and a particular geometri- 
cal setting capable of eliminating the NXLD signal and 
leaving in evidence the nonreciprocal dichroic effect. 

Even accepting the fact that by a quite remote chance 
the authors in Ref. 15 had fallen on the correct setting, 
there still remain some paradoxical aspects with their 
result. The group PI' ja is ME; however a recent mea- 
surement of the ME effect 19 has given a negative result. 
Moreover this group does not seem to be among the pos- 
sible candidates of the magnetic symmetry of the ground 
state in the AFI phase, as obtained by the minimiza- 
tion of an effective Hubbard Hamiltonian with degener- 
ate bands. 11-13 Finally, the reduction to the group P2'/a 
seems not compatible with the features of the Bragg- 
forbidden reflection (lll) m : the signal is present only at 
the quadrupole energies (5465 eV) and it is absent at 
the dipole energies (5470-5490 eV), 9 while no extinction 

rules are present at the dipolar energies for the P2' /a 

39 

group. 

Faced to these paradoxes, we have tried to pro- 
vide an alternative interpretation of the "nonreciprocal" 
results, 15 invoking a geometrical origin for the reversal 
of the linear dichroism with the external magnetic field. 
However the question can only be settled by experimen- 
tal tests in which one has the possibility to check the 
twin after each transition to the AFI phase, by means 
of an x-ray scattering equipment: this would give a full 
geometrical control of the system to be analyzed. Given 
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the implications of the nonrcciprocal effect, if confirmed, 
for the physics of V2O3, its unambiguous experimental 
determination is of the outmost importance. 



We would like to acknowledge interesting discussions 
with L. Paolasini on his measurements of the Tcmpleton 
scattering in V2O3. 
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